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We study ultracold atoms in an optical lattice with two local minima per unit cell and show that 
the low energy states of a multi-band Bose-Hubbard (BH) Hamiltonian with only pair-wise inter¬ 
actions is equivalent to an effective single-band Hamiltonian with strong three-body interactions. 

We focus on a double-well optical lattice with a symmetric double well along the x axis and single 
well structure along the perpendicular directions. Tunneling and two-body interaction energies are 
obtained from an exact band-structure calculation and numerically-constructed Wannier functions 
in order to construct a BH Hamiltonian spanning the lowest two bands. Our effective Hamiltonian 
is constructed from the ground state of the iV-atom Hamiltonian for each unit cell obtained within 
the subspace spanned by the Wannier functions of two lowest bands. The model includes hopping 
between ground states of neighboring unit cells. We show that such an effective Hamiltonian has 
strong three-body interactions that can be easily tuned by changing the lattice parameters. Finally, 
relying on numerical mean-field simulations, we show that the effective Hamiltonian is an excellent 
approximation of the two-band BH Hamiltonian over a wide range of lattice parameters, both in 
the superfluid and Mott insulator regions. 

PACS numbers: 67.85.-d, 37.10.Jk, 03.75.Lm 


I. INTRODUCTION 

Experiments with ultracold atoms in optical lattices 
have been used to realize the Bose-Hubbard (BH) model 
[I]. Its interaction driven quantum phase transition was 
first observed in §2]. Ultracold atoms in optical lat¬ 
tices are highly tunable systems and are now increas¬ 
ingly used to simulate other quantum many-body Hamil¬ 
tonians M- In recent years, Hamiltonians with com¬ 
plex band structure have been experimentally realized, 
such as double-well lattices [SHE], honeycomb, triangu¬ 
lar, one-dimensional stripe and Kagorne lattices mm, 
and artificial graphene HJ. These experiments have ex¬ 
tended the scope of the BH model to include excited 
bands and richer on-site interactions. Such Hamiltoni¬ 
ans are being used to study collective phenomena ranging 
from modified quantum phases mm3, topological mat¬ 
ter 0 m ns ns, and to the formation of long-lived un¬ 
conventional Bose-Einstein condensates in excited bands, 
which were theoretically studied in Refs. mm 3 and ob¬ 
served by Refs [2D [22] ■ 

Recent studies have shown that couplings between ad¬ 
jacent unit cells due to atom-atom interactions are also 
important and can lead to density-induced tunneling. For 
cubic lattices, such additional terms in the BH Hamilto¬ 
nian have been shown to measurably modify the location 
of the superfluid to Mott phase boundaries [25] [23 • 

Many-body Hamiltonians with three-body interaction 
are also being explored. Three-body interactions can lead 
to significantly modified quantum phases [251127] . and 
can yield Pfaffian-like ground states [25H5T] that have 


promising uses in quantum computation , 521 - Following 
the proposals by Refs. [TTTTHTT(Tj . the first experiments con¬ 
firmed the presence of weak effective three-body inter¬ 
actions in optical lattices [3714401 . These effective three- 
body interactions, due to virtual transition of atoms to 
higher bands, are smaller than the usual two-body inter¬ 
action. Recent proposals have investigated ways to re¬ 
duce the strength of this two-body term, thereby enhanc¬ 
ing the role of three-body interactions. In fact, Ref. @13 
by driving the lattice at rf frequencies and Ref. @2] by 
adding resonant radiation to couple internal states of an 
atom have proposed ways to turn off the two-body in¬ 
teraction altogether. It was shown in @3] that arrays 
of superconducting Josephson junctions can also mimic 
strong three-body interaction. 

In this paper we show that very strong effective three- 
body interactions can be created using optical lattice po¬ 
tentials with two local minima per unit cell. We study 
trapped atoms in a double-well optical lattice with three 
or less atoms per site, and show that over a wide range 
of lattice parameters, the on-site interactions can be de¬ 
scribed by an effective Hamiltonian of the form 

3 1 

H eS = V —r (l) 

ml 

m=1 

where b^(b) are creation(annihilation) operators, and F m 
represents the effective m-body interaction energy. We 
show that tuning the lattice parameters can create situa¬ 
tions when r 3 ^>r 2 ,i.e.,a system with a strong effective 
three-body interaction energy. 

The remainder of the paper is setup as follows. In 
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FIG. 1: (color online) (a) Contour plot of the optical lattice 
potential in the xy plane, where the potential minima are 
in dark blue. The white box encloses a unit cell of length 
a and a/2 along x and y, respectively. Each unit cell has a 
double well along the x axis, labeled L and R, and a single 
well along the y and 2 axes, (b) Schematic of a symmetric 
double-well potential (thick solid black curve) in a unit cell 
along the x axis. The barrier height between the left and 
right wells is B. The Wannier functions for the ground and 
first-excited band, wi(x) and W 2 {x), are shown as blue curves. 
They are computed for a lattice with Vq/Er = 17.0, Vi/Vo = 
1.3, V 2 /Er = 40.0 and fcfcn = 7t/4. The two horizontal red 
curves represent their energies with splitting 5. 


Sec. |TT] we introduce the double-well optical lattice po¬ 
tential and numerically obtain the band structure and 
Wannier functions. We then introduce the multi-band 
BH Hamiltonian and calculate hopping and interactions 
parameters of this Hamiltonian. In Sec. |III| we per¬ 
form a mean-field calculation to determine its phase di¬ 
agram and discuss the Superfluid-Mott insulator tran¬ 
sition. The effective Hamiltonian model is introduced 
in Sec. IV and constructed in a three-step process. In 
Sec. |IV A we discuss the initial step of obtaining a many- 
particle basis and energy levels for each unit cell. The 
effective Hamiltonian in a unit cell of the form ([I]) is de¬ 
termined in Sec. IV B| We also show how Ts/^ can be 
optimized by tuning the lattice parameters. In Sec. |IV C| 
we study the coupling between the many-particle states 
in adjacent unit cells and introduce the effective tunnel¬ 
ing Hamiltonian. We discuss the validity of the effective 
Hamiltonian by comparing the numerical mean-field re¬ 
sults obtained for the effective Hamiltonian picture and 
the full BH Hamiltonian in Sec.[Vj Finally, we summarize 
our results and present conclusions in Sec. |Vl| 


II. HAMILTONIAN FOR THE DOUBLE-WELL 
OPTICAL LATTICE 

We are interested in atoms in optical lattice potentials 
that have two nearly degenerate local minima in an unit 
cell. In particular, we focus on a lattice with a double¬ 
well structure along the x axis and single-well structure 
along the perpendicular y and 2 axes. Such a lattice can 



FIG. 2: (color online) Numerical parameters for the BH 
Hamiltonian for 8 ‘ Rb on a log-linear scale (a) and linear-linear 
scale (b) plotted as functions of lattice depth Vo in units of 
Er , at V 1 /V 0 = 1.3, V 2 /ER = 40.0 and bk,L = tt/ 4. The tun¬ 
neling energies Ji, —J 2 and J± are shown as dot-dashed blue 
curves. The band-gap S is plotted as a solid red curve. The 
interaction energies U r = U\, U 2 and U 12 are shown as dashed 
black curves. They are calculated for a scattering length of 
a 3 = 100ao, where ao = 0.0529 nm is the Bohr radius. 


be constructed by using a laser with wave vector kL and 
its first higher harmonic. We focus on 

V(x) = —Vo cos 2 {kLx) — Vi cos 2 [2 Ul{x + &)] 

- V 2 [cos 2 (2 k L y) + cos 2 (2fc L 2)] , (2) 

where Vo 1 2 are lattice depths and & is a parameter 
that determines whether the lattice has a symmetric or 
asymmetric double-well. The lattice has a periodicity of 
a = ir/kL along x axis, and a/2 along the perpendicu¬ 
lar directions. Throughout, we express energies in units 
of the recoil energy along the x axis, Er = h 2 kf j /2m a , 
where m a is the atomic mass. We note that the recoil en¬ 
ergy along the perpendicular directions is Er± = 4 Er. 
Figure [Ta| shows a contour plot of the optical-lattice po¬ 
tential in the xy plane, while Fig. |lb| shows a symmetric 
double-well for = 7r/4 along the x axis. 

In this paper, we solely focus on the symmetric double¬ 
well lattice with k^b = 7r/4 and allow the lattice depths 
to vary. The barrier height between the left and right 
wells is B = (Vi — V(j/4) 2 /Vj. The barrier disappears 
when Vi/Vo <1/4. In the limit Vi/Vo —¥ 00 , the 
potential approaches a single-well potential with period 
a/2 along all directions. Some remarks about the case 
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&l 6 ^ 7 r /4 are made at the end of Sec. IV B 


The band structure and corresponding Bloch func¬ 
tions are obtained numerically. They are independently 
computed for the x , y and z directions and lead to a 
band dispersion e s (k) = e XtCt:D (k x ) + ey, ay (ky) + e Zta *(k z ), 
where a = ( a x , a y , a z ) are the band indices and k is the 
quasi-momentum in the first Brillouin zone. The low¬ 
est two bands, a x = 1 , 2 , along the x axis are quasi¬ 
degenerate for deep lattices. Along the perpendicular 
directions, the ground band a y = a z = 1 is far removed 
from the first excited band. The Wannier functions, la¬ 
beled by unit cell i = (i x ,i y ,i z ) and a, are a product 
of one dimensional Wannier functions, one for each axis, 
i.e., w it s(x) = wi ax {x)w hay (y)wi az (z). These functions 
can be obtained with the Marzari-Vanderbilt scheme [44] 
of constructing maximally localized Wannier functions. 
This was used for a double-well optical lattice in Ref. [55]. 
Here, we follow Refs. mm and Wannier functions are 
found as eigen states of the position operators x, y and 
z within the subspace of all the Bloch functions with 
band index a. To ensure real-valued Wannier functions, 
a real-valued discrete variable representation with peri¬ 
odic boundary conditions m over M = (M x , M y , M z ) 
unit cells is used. Figure [Tb] shows the Wannier func¬ 
tions Wi,a x {x) for the lowest two bands along the x axis 
for a symmetric lattice. We see that the Wannier func¬ 
tions are exponentially suppressed between unit cells, but 
that they are extended over the barrier between the left 
and right well. This is typical for the lattice parameters 
that are the focus of this paper. 

Since the lowest two bands along x axis are quasi¬ 
degenerate for deep lattices, we construct a multi-band 
Bose-Hubbard (BH) Hamiltonian H spanning the low¬ 
est two bands along the x axis, and only the ground 
band along the perpendicular directions. For simplic¬ 
ity, we drop the band indices a y = a z = 1 , and use 
a = 1,2 to refer to the band index along the x axis. 


We find H = 


II- 


cell 

neighbor hopping Hamiltonian is 


Ei {#i hop 


, where the nearest- 


jjhop __ _ 


Ji- ^ ' (^o^i+luO f f ^ C ’J 

c*ei,2 

'y ] J a ^ a i. Q a i+l„..a + h.C.^j . (3) 


ael ,2 


The operators a\ a and a\ M create and annihilate a par¬ 
ticle in the Wannier function of unit cell i and band 
a = 1,2, respectively. The abbreviation h.c. is the 
hermitian conjugate. The location i + l x denotes the 
unit cell (i x + 1, i y , i z ). A similar notation is used for 
the other directions. The nearest-neighbor (NN) tun¬ 
neling energy in band a along the x axis is J a = 
(1 /M x )J2 k cos {ka)e x , a (k x ), where the sum is over the 
M x allowed quasi-momenta k x . Along the perpendicular 
directions the NN tunneling energy in the ground band is 
J±. There is no cross-tunneling between separate bands 
as the Wannier functions are orthogonal. We have not 


included next-nearest neighbor tunneling since our nu¬ 
merical computation show that they are nearly two or¬ 
ders of magnitude smaller than the NN tunneling for the 
lattice depths on which we focus in this paper. 

The Hamiltonian ffP 11 for unit cell i is 



+ 2 t/ 

+ ^U 12 (4,i a i,i a i,2 a i,2 + 4,2 a i,2 a i>i a i,i) > ( 4 ) 

where hi iQ , = a\ a a i iQ and the band-gap between bands 
a = 1 and 2 is 6 = (1 /M x )J2 kx { € x,2(k x ) - e xA (k x )}. 
The on-site interaction terms U\, U 2 and U\ 2 are ob¬ 
tained numerically from the local Wannier functions as 
U a = g f Wi,a(x)wi !a (x)wi t a(x)wi t a(x)dx and U a p = 
g f w a (x)w a (x)wp(x)wp(x)dx, where g = 47 t h 2 a s /m a 
and a s is the s-wave scattering length. As a consequence 
of the nature of the Wannier functions, the term U\ 2 is 
comparable to U\ and U 2 giving rise to strong density- 
density and pair-tunneling terms. It is the interplay be¬ 
tween this pair-tunneling term and the band gap contri¬ 
bution that will enable us to achive large effective three- 
body interaction. 

For a symmetric double-well lattice, other atom-atom 
interaction terms between Wannier functions in bands 1 
and 2 within a unit cell are strictly zero due to the parity 
of the Wannier functions. For small asymmetries, these 
other terms are small and can be neglected. Finally, we 
have not included terms due to atom-atom interactions 
between Wannier functions in neighboring unit cells, be¬ 
cause they are an order of magnitude or more smaller 
than those in Eq. J4|. 

Figures [2a| and [2*b| show parameters of the BH Hamil¬ 
tonian as a function of lattice depth Vq for 87 Rb atoms at 
a fixed ratio of V\ / Vq and fixed V 2 . The plot shows that 
the interaction energies U\, U 2 and U\ 2 are of equal im¬ 
portance, satisfy U 2 > XJ\ 2 > XJ\, and are of the order of 
the band gap 5. The tunneling energies are much smaller 
than the interaction energies. The tunneling energy J 2 
for band 2, as expected, is negative and \J 2 \ > \ J\\. Our 
choice of V 2 /Er = 40.0 is such that, as shown in the 
next section, our Hubbard Hamiltonian has a superfluid 
region. 


III. PHASE DIAGRAM USING DECOUPLING 
APPROXIMATION 

In this section, we obtain the mean-field phase diagram 
for the BH Hamiltonian. A knowledge of the quantum 
phases as a function of lattice parameters is essential to 
prepare the system in the required many-body ground 
state. In addition, the results will be used in Sec.[V]to dis¬ 
cuss the validity of the effective Hamiltonian model. The 
phase diagram is obtained in the mean-field decoupling 



4 



FIG. 3: (color online) Phase diagram for 8 ' Rb in a symmetric 
double-well lattice as a function of lattice depth Vo in units 
of E r at V 2 /E r = 80.0, Vi/Vo = 1.3 and scattering length 
a a = 100ao. The dark blue lobes are the Mott lobes for 
N = 1, 2 and 3 atoms per unit cell. All other regions represent 
SF phases, with SI denoting the case 0 ^ 0, 0 = 0 and S 2 
denoting ■02 0, 0 =0. 

approximation |35j by introducing a real-valued homoge¬ 
neous superfluid order parameter for each band a = 1,2, 
i.e., ip a = (a\ a ) = (a^o) and approximate a\ a a jjQ by 

ip a (°i Q + Oj,a) — V’a for i ^ j- In the grand-canonical 
ensemble, this leads to the grand-canonical potential or 
Hamiltonian H MF = i MF — H a \ a a i,a), where 

HfF =i?r n - (a\ a + a 1>a ) (2 J a + Z±J±) 10 

aG 1,2 

+ (2 J a + z±J±) ipa (5) 

ae 1,2 

and £_l is the number of neighboring unit cells along the 
perpendicular directions. In our case, z± = 4. The 
ground state of the grand-canonical potential H MF at 
a given chemical potential n corresponds to order pa¬ 
rameters (0i, 02 ) that minimize its ground state energy. 
We have performed these calculations in a basis of Fock 
states | rii, 77 , 2 ), where n a is the number of atoms in band 
a in an unit cell. We reach numerical convergence for 
ni,ri 2 < -/V max = 10. A Mott insulator (MI) state corre¬ 
sponds to 0 = 0 = 0, while a superfluid (SF) state is 
obtained when either 0i, 02 7 ^ 0. 

Figure [3] shows the mean-field phase diagram for our 
symmetric double-well lattice. The figure shows Mott 
lobes up to N = 3 atoms per unit cell. We find that the 
Mott lobes with odd N are reduced in size compared to 
adjacent lobes with even TV, a feature also noted and dis¬ 
cussed in D3I- For the SF phase 51, the order parameter 
ipi 0, but tp 2 = 0, while for the SF phase 52 0 7 ^ 0 
but 0 = 0. Both the order parameters have a jump at 
the 51-52 phase-boundary corresponding to a first-order 
phase transition. For the relatively shallow depths used 
in all directions in Fig. [3j superfluidity extends over all 
dimensions. For larger depths along either x or y and 
z, the superfluid can be localized to fewer dimensions. 


A further stability analysis around the mean-field solu¬ 
tion is then required to identify this distinction j49ti51j . 
Such knowledge, however, is not crucial for deriving the 
effective Hamiltonian model. 

A cutting off of the Mott lobes for N = 1 and 2 at the 
phase-boundary to phase 52 is apparent in Fig. [3| From 
a separate analytical calculation of the phase-boundary 
of H mf , excluding the interaction induced pairing terms, 
we find that this cut-off occurs when 20 + z±J± =0. 
For the lattice parameters used in Fig. [3]this occurs when 
Vo/ Er ss 20, consistent with our numerical simulations 
including the pairing term. 

The lattice parameters in Fig. [3] were particularly cho¬ 
sen so as to have a good representation of all the three 
quantum phases. Since our primary aim in this paper is 
to show the presence of large effective three-body inter¬ 
actions in double-well optical lattices, we do not provide 
phase diagrams for a wider set of lattice parameters. In¬ 
stead, we note that in our case, we are mainly interested 
in the phase 51 and the corresponding SF-Mott phase¬ 
boundary. Also, the region 51 grows relative to those of 
other phases as we decrease the lattice depth along the 
perpendicular directions. 

IV. EFFECTIVE HAMILTONIAN 

We now show that the low energy states of the multi¬ 
band BH Hamiltonian defined in Eqs. ^ and © can 
be obtained from an effective Hamiltonian H e ff that has 
strong three-body interactions. Constructing H c g is a 
three-step process. The first step is discussed in Sec. |IV A| 
and involves diagonalizing the on-site Hamiltonian _f/j Cen 
in unit cell i to obtain many-particle (MP) energy lev¬ 
els. Using these MP levels an effective on-site interac¬ 
tion Hamiltonian is constructed in Sec. |IVB| Finally, in 
Sec. |IV C| we calculate the effective tunneling Hamilto¬ 
nian that couples the MP states of adjacent unit cells. 

A. Many-particle energy levels 

We are interested in lattice parameters for which the 
average band gap S and interaction energies U\,U 2 and 
t/12 are much larger than the tunneling energies J \, J2 
and J±. We thus diagonalize the on-site Hamiltonian 
7J. ce11 in the Fock state basis |ni, N — n\) for N atoms in 
unit cell i with ni atoms in band a = 1. This gives the 
many-particle (MP) eigen energies £^ and eigen states 
\v, N) with v = {1,..., iV + 1}. We refer to these v states 
as vibrational states. The MP eigen-states in terms of the 
Fock states are 

N 

\v,N)=Y l C%>mni,N-n 1 ), ( 6 ) 

ni—0 

with real coefficients C^i\N). It is good to note that the 
many-body ground state of the Hamiltonian in Eq. ([ 3 ]) 
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FIG. 4: (color online) (a) Many-particle energy levels £^ 
for two atoms per unit cell, as functions of lattice depth Vo in 
units of Er. The solid blue lines are for Vi/Vo = 1.3, while the 
dot-dashed magenta curves are for V\/Vq = 1.7. (b) Similar 
plot for three atoms per unit cell. Plots are for s 'Rb in a 
symmetric lattice with V 2 /Er = 40.0 and scattering length 
a s = 100 ao. 


in a Mott lobe is a product of | is,N) states with v = 1, 
one for each site. 

For N = 1 the eigen states are simply the Fock states 
11,0) and |0,1) with energy £^ = —5/2 and £^ = 5/2. 
A plot of the band gap 5 as a function of lattice depth 
Vo is shown in Fig. [2a] Figures [4a] and [4b] show the 
MP energy levels £^ cts a function of lattice depth Vo 
and at fixed Vi/V 0 and V 2 for two and three atoms, re¬ 
spectively. Their behavior can be understood by not¬ 
ing that the band gap <5 decreases with increasing lat¬ 
tice depth Vo, and for larger depths 5/U 1 with 

U « U\ ~ U 2 ~ U 12 - Consequently, we can show using 
perturbation theory that for N = 2, E^ = —5 2 /(2U), 


E^ = 2U and £ 2 
7(1) 


(3) 


= 2U + 5 2 j(2U). For N = 3 atoms, 


. 3 =2U - 5 - 35 2 /6U, £f> = 2U + 5 - 35 2 /6U, and 


2 b 


£3 = £3 = 6C7 + 35 2 /6U. In conjunction with Fig. ^ 

these expressions give a fair estimate of £^ for large Vo. 
We also note that the energy of the ground-vibrational 
v = 1 state is well separated from that of the state v = 2 
with a spacing much larger than Ji, J 2 or J±. For the 
effective Hamiltonian picture, we assume that the atoms 
only populate the v = 1 states. 



FIG. 5: (color online) Log-linear plot of the ratio of three- to 
two-body interaction strength r 3 /r 2 as a function of lattice 
depth Vo in units of Er for various ratios of V] /Vo. The 
horizontal red line denotes r 3 /T 2 = 1. The plot is for 8, Rb 
in a symmetric lattice at Vi/Er =40.0 and scattering length 
o s = IOOoq. 


We should further note that interactions in Eq. Q mix 
the Fock states \ni,N — n\) and the vibrational ground 
states \u = 1, N) at large lattice depths have considerable 
contribution from Fock states with atom population in 
the first excited band. 


B. Effective interaction Hamiltonian 

We can now construct the effective on-site interaction 
Hamiltonian similar to Eq. ([!]) in unit cell i based on the 
large energy separation between the v = 1 and 2 vibra¬ 
tional states for each N. In other words we assume that 
the atoms only populate the ground vibrational states 
with energies £^ that are reproduced by the effective 
on-site interaction Hamiltonian 

fi eff=EEi r 4T, (7) 

lit. 

1 m= 1 

where bi and b\ are effective bosonic annihilation and 
creation operators for unit cell i and state v = 1. The 
coefficients T m are the m— body interaction strengths, 
which we restrict to m < 3, and are found by mapping 
£^ to the eigen-energies of . Some algebra shows 
T 1 = e[ 1} , r 2 = 4 1) -2£ 1 (1) and r 3 = £^ 1} -3(£^ } -£ 1 (1) ). 

Figure [5] shows r 3 /r 2 as a function of lattice depth Vo 
for various V\/Vq. We observe that r 3 /r 2 increases with 
Vo and V 1 /V 0 , i.e. with increasing lattice depth. In fact, 
choosing appropriate lattice parameters can produce a 
three-body strength that is larger than the two-body one. 
Using the same approximations as used in Sec. |IV A| for 
large lattice depths, we find r 2 «U 5 2 /{2U) + 0(5 3 ) 
and r 3 « 2U — 2.5(5 +1.3d 2 + (!I((5 3 ). Consequently, 5 and 
r 2 —> 0 while r 3 remains finite for larger lattice depths 
and r 3 » r 2 . 

For an asymmetric double-well lattice with k^b 7 ^ 7 t/ 4 , 
5 and equivalently T 2 depend on the asymmetry and re- 
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FIG. 6: (a) Plot of ratio 1— Anm / y/(N + 1 )M as a function 
of lattice depth Vo. The coefficient Aoi is not shown because 
the corresponding ratio is zero, (b) Plot of B pq as a function 
of lattice depth Vo. Both the panels are for s 'Rb with k^b = 
7t/4, Vi/Vo = 1.3, Vi/Er = 40.0 and a s = 100ao. 


main finite for a large lattice depth. A symmetric lattice 
is therefore the most favorable case with which to reach 
sizeable T 3 . We should also note that it is not possible 
to obtain r 2 <C T 3 for a single-well lattice since the band 
gap increases with lattice depth. 


C. Effective tunneling Hamiltonian 

In this section we study the coupling between the MP 
states u, N) of neighboring unit cells. Since tunnel¬ 
ing involves two adjacent unit cells, we label states as 
\v, N) i|z/, M) where N and M are the atom numbers 
in unit cell i and i -P1^ , i + ly or i + l z , respectively. As in 
the previous subsection, we restrict our discussion to the 
v = 1 ground state and simply write the initial and final 
states as \N, M; i, i + 1}} = \v = 1, N)\\v' = 1, 

With these definitions we can express tunneling be¬ 
tween neighboring unit cells along the x axis as 

i N,M 

IAf+1, M— 1; i, i + 1 X )){{N , M; i, i + l^j + h.c. j- 
where the sums over N and M are restricted to 0 < N < 


M , and using Eqs. ([s]) and 


-JiAnm = «iV-|-l,M—l;i, i + l x |i? i hop |A', M; i, i + l x )} 

N M 

\Z(n + l)mCp(N) 

^n+1 (N + 1) 


n—0 m—1 


X C«(M)C« 1 (M + 1), 


(9) 


are unit-cell-independent dimensionless real coefficients. 
Here, we only include tunneling through the term propor¬ 
tional to J\ in Eq. ([3]), since the v = 1 vibrational level 
is predominantly determined by the \N, 0) Fock state. 
Thus, the contribution from tunneling through the J 2 
term is negligible. We limit ourselves to N,M < 3, i.e. 
up to three atoms per site, and only require the six in¬ 
dependent Aoi, A 02 , A 12 , A 03 , Ais and A 2.3 ■ In Fig. [6a| 
we plot the ratio 1 — Anm/^J(N + 1 )M for these six 
coefficients as a funct ion of lattice depth Vo- The ex¬ 
pression yj(N + 1 )M corresponds to the coefficient ex¬ 
pected for a simple single-well BH Hamiltonian. The 
coefficients Anm deviate only slightly from these values 
for Vo/En < 30 in Fig. 6a More generally, we find that 
the correction becomes larger for larger lattice depth. 

The effective tunneling Hamiltonian has atom- 

number dependent parameters. We find it convenient to 
Wrfte H x°eS in an alternative form as a sum of multi-body 
hopping operators 


rrhop _ _ 


k «-l 


A HEM 6 !) 0-u) 

i Q,P 

x (6i) p +/i.c.j, (10) 

where the sums over p and q are restricted to 0 < p < q. 
An ?n-body hopping operator corresponds to terms where 
p+q = m. The six coefficients B pq are found from solving 
the linear system 


3 <?—1 

Anm = EE 8 « E{N,M,p , q) 

q—1 p —0 


(ii) 


With F(N, M,p,q) = y/(N+l)M ^J M q) 


Figure 6b shows a plot of B pq as a function of lat¬ 


tice depth Vo- The most important term is Bqi = 1, 
the usual single-particle tunneling term between adja¬ 
cent unit cells. The next two leading terms are the two- 
and three-body operators with strength B 02 and B i 2 , re¬ 
spectively. Their values are negative. Other terms are a 
factor of two or more smaller and positive. 

Figure [4b] shows that the energy separation between 
the v = 1 and 2 vibrational states for N = 3 atoms de¬ 
creases with increasing ratio V\/Vq. This might seem to 
invalidate the assumption that the atoms only populate 
the ground vibrational state for larger lattice depths. For 
such a scenario with large V\/Vo, however, we can choose 
a smaller lattice depth Vo to ensure that the ground vi¬ 
brational state is well separated, and still have a large 
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assuming that the atoms only occupy the ground vibra¬ 
tional (v = 1 ) states and coupling to excited v states 
is negligible. This assumption is justified as the energy 
gap between the ground and excited vibrational states 
(see Figs. 4a and 4b) is much larger than the coupling 
strengths, which are of the order of the tunneling ener¬ 
gies. We now verify this assumption further by compar¬ 
ing mean-field phase diagrams. 


For our purposes it is sufficient to further simplify H e g- 
to 



Lattice Depth Vq/Er 


FIG. 7: (color online) (a) Plot of the mean-field ground-state 
energy E g as a function of lattice depth Vo in units of Er and 
V 2 / Eit = 40 for different chemical potentials p. The solid 
lines represent the results for the full Hamiltonian, while the 
dashed lines show those for the effective Hamiltonian. Lat¬ 
tice parameters are chosen such that the system is in the SF 
region SI. (b) Mean-field SF-Mott phase diagram for H c g as 
a function of lattice depth Vo/ Er at V 2 / Er = 80. The deep 
blue areas are Mott lobes surrounded by a SF region. Other 
parameters are the same as those used in Fig. [3] 


F 3 > r 2 . Moreover, the coupling between the ground 
and excited vibrational states in neighboring unit cells 
are of the order of the tunneling energies with strength 
proportional to Jf/y/U 2 + S 2 , which decrease with in¬ 
creasing V\ /Vo, and stay much smaller than the energy 
separation between the ground and excited vibrational 
states. 

The effective tunneling Hamiltonians along perpendic¬ 
ular directions are obtained by replacing J\ by J±. The 
numerical values for the coefficients B pq are the same as 
those along the x axis. This finally leads to the effective 
tunneling Hamiltonian H^ p = H^° p s +Hy° p s +H^° p s and 

the total effective Hamiltonian H e s = + H^g P . 


V. VALIDITY OF THE EFFECTIVE 
HAMILTONIAN 


In this section, we discuss the validity of the effective 
Hamiltonian developed in Sec. [TV] It was obtained by 


i ^ m =1 

-J± (b\bi + i y +^ 6 i + i^ 


Jib\b i+1:D 
+ ft.c.j , 


( 12 ) 


where we have only kept the leading-order single-body 
tunneling terms in Eq. (10). Thus, H e s is a single¬ 
band BH Hamiltonian with an added effective three-body 
term proportional to T 3 . Similar to Sec. |III| we intro¬ 
duce a real order parameter ip = (b i) = (&() and deter¬ 
mine the ground-state energy, E g , for the corresponding 
mean-field Hamiltonian in the grand-canonical ensemble 
by self-consistently minimizing the energy with respect 
to the order parameter ip. 


We first confirm that H e g is a good model for lattice 
parameters where the full Hamiltonian in Eqs. (J3| and 
Q has a superfluid ground state in the 51 region. Fig¬ 


ure 


7a compares the ground-state energy E g of the full 


mean-field and the effective Hamiltonian as a function of 
lattice depth Vo for various values of p. This graph is 
computed for V^/Er = 40, half the size used in Fig. [ 3 ] 
The mean-field ground state of the full Hamiltonian is 
then in the 51 region, with order parameter ip 2 = 0, 


for the entire parameter space (/it, Vo) shown in Fig. 7a 


Similarly, the ground state of the effective Hamiltonian 
is superfluid. The mean atom number per unit cell for 
the four /is shown in Fig. [7a]are n = 1.56, 2.37, 3.19 and 
4.01 at a lattice depth of Vq/Er = 30, respectively. 


The ground-state energies for the two models show 
similar trends, have differences that are nearly indepen¬ 
dent of Vo and agree better for smaller chemical poten¬ 
tials. The larger discrepancy for larger n is not due to 
virtual excitations to higher v vibrational states or m > 1 
hopping terms since J\B pq or J\B pq with p, q ^ 0,1 are 
much smaller than the recoil energy. For n > 3 we need 
to include four- or higher-body interaction terms to get 
better agreement. 


Next, we study whether H e g is valid as we approach 
the SF-Mott phase-boundary. Unlike the ground state 
energy, the mean-field phase boundary can be analyti¬ 
cally obtained within second-order perturbation theory 

[48] . treating — (2Ji + z± J±) (b\ +b-A ip as a perturba- 
















tion, and is given by 


n t fi i_ (g + 1) _. 

eff 1 + -r 1 + ig(r 3 - 2 t 2 ) - y*r 3 + 

g'A>:7 1 

(r 3 - r 2 + ro + y(2T 2 - 3r 3 ) + i<? 1 2 3 4 5 6 7 8 9 10 11 12 r 3 - ? J ’ 

with J e ff = (2Ji + z± J±) and the Fock state |g) is the 
zeroth-order ground state of the grand-canonical poten¬ 
tial. Figure [7b| shows the corresponding SF-Mott phase 
diagram for V 2 /E R = 80 and can thus be directly com¬ 
pared with the phase diagram shown in Fig. [3] for the 
full Hamiltonian. For Vq/E r > 20 the Mott lobes for 
N = 0, 1 and 2 are nearly identical. The N = 3 lobe for 
the effective Hamiltonian, however, is significantly larger 
than that of the full Hamiltonian. In fact, the tip of the 
Mott lobe has shifted to smaller lattice depth. We have to 
include T 4 to correctly model this lobe. Finally, our effec¬ 
tive Hamiltonian does not describe S2 for Vq/E r < 20. 
In this phase, the full Hamiltonian has a non-zero order 
parameter in the second band. Our effective Hamiltonian 
assumes no population in the second band and can not 
represent these states. 

VI. CONCLUSION 

We have shown that the low energy states of a sys¬ 


tem of trapped atoms in a double-well optical lattice can 
emulate a Hubbard model with strong three-body inter¬ 
actions. The optical lattice with a double-well potential 
along one direction and single period lattice along the 
remaining directions has two nearly degenerate bands. 
The Hamiltonian has a strong pair-tunneling contribu¬ 
tion between the two Wannier functions within a unit 
cell. The interplay between this interaction and the band 
gap plays an important role in determining the behavior 
of the system. In particular, we find that the low energy 
states of such a system are quite accurately described by 
an effective single-band Hamiltonian with a strong three- 
body interaction, whose strength can be easily tuned by 
changing the lattice parameters. Surprisingly, tunneling 
between the neighboring unit cells in the effective Hamil¬ 
tonian has, to good approximation, the same structure as 
that for a particle hopping in a single-band BH model. 
By comparing the numerically obtained ground state en¬ 
ergy and phase diagram of the full and effective Hamil¬ 
tonians, we verified that the effective Hamiltonian model 
is an excellent approximation over a wide range of lattice 
parameters, both in the superfluid and Mott insulator 
phases. 
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